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Numerical Modeling of Pulse Detonation Rocket 
Engine Gasdynamics and Performance 


C. I. Morris* 

NASA/George C. Marshall Space Flight Center 
TTMO/Propulsion Research Center 
Marshall Space Flight Center, AL USA 35812 

Pulse detonation rocket engines (PDREs) offer potential performance improvements over conventional 
designs, but represent a c h all eng i n g modeling task. A quasi- 1-D, finite-rate chemistry computational fluid dy- 
namics model for PDREs Is described and implemented. Four different PDRE geometries are evaluated in this 
work: a baseline detonation tube, a detonation tube with a straight extension, and a detonation tube with two 
types of converging-diverging nozzles. The effect of extension length, and C-D nozzle area ratio, on the gasdy- 
namics and performance of a PDRE is studied over a wide range of blowdown pressure ratios (1-1000). The 
overall effect of straight extensions and converging diverging nozzles on performance and blowdown time is 
compared. The results are also compared to a steady-state rocket system using similar modeling assumptions. 


Introduction 

Pulse detonation engines (PDEs) have generated con- 
siderable research interest in recent years as a chemical 
propulsion system potentially offering improved perfor- 
mance and reduced complexity compared to conventional 
gas turbines and rocket engines. The detonative mode of 
combustion employed by these devices offers a theoretical 
thermodynamic advantage over the constant-pressure de- 
flagrative combustion mode used in conventional engines. 
However, the unsteady blowdown process intrinsic to all 
pulse detonation devices has made realistic estimates of the 
actual propulsive performance of PDEs problematic. The 
recent review article by Kailasanath 1 highlights some of 
the progress that has been made in comparing the available 
experimental measurements with analytical and numerical 
models. 

The pulse detonation rocket engine (PDRE) 2 has re- 
ceived a comparatively more limited research interest than 
the pulse detonation engines intended for airbreathing use. 
While conceptually similar, the rocket application of this 
technology requires investigation of the performance char- 
acteristics over a wide range of ambient pressures. In con- 
ventional rocket engine designs, the converging-diverging 
(C-D) nozzle plays a critical role in converting the ther- 
mal energy of the combustion products in the combus- 
tion chamber to directed kinetic energy in the exhaust. 3 
Thus, there has been considerable recent research interest 
in studying the use of nozzles for enhancing PDE perfor- 
mance. 4 No attempt is made here to provide an exhaustive 
review of nozzle research for PDEs, but the most recent 
results, relevant to the rocket application and the current 
study, should be listed. Cooper and Shepherd 5 have re- 
cently developed a partial-fill model for modeling the effect 
of straight extensions and diverging nozzles on detonation 
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tubes at atmospheric conditions. This model makes use 
of the available experimental impulse data in the litera- 
ture, and the partial-fill numerical simulations of Li and 
Kailasanath. 6 Numerical studies of the effects of nozzles 
on a single detonation tube over a wider range of ambient 
pressures have been performed by Cambier 7 and Mohanraj 
and Merkle. 8 A study incorporating multiple tubes into a 
common nozzle was performed by Mohanraj et al. 9 Tal- 
ley and Coy 10 have developed a 0-D, constant- y model for 
studying the performance of an idealized constant- volume 
combustor with a nozzle over a wide range of pressure ra- 
tios. 

In a previous paper by the author, 11 parametric studies 
of the single-shot performance of a single, straight deto- 
nation tube were reported over a wide range of pressure 
ratios. A 1-D, unsteady method of characteristics code, 
employing a constant- y assumption behind the detonation 
front, was developed for that study. Models of this type 
are computationally inexpensive, and may be used to illus- 
trate the critical role of nozzles in the rocket PDE appli- 
cation. For example, a plot showing the specific impulse 
of various PDRE and steady-state rocket engine (SSRE) 
configurations as a function of blowdown pressure ratio is 
shown in Figure 1. The blowdown pressure ratio is de- 
fined throughout this work as the ratio of the initial fill (or 
combustion chamber) pressure, P c j , to the ambient pres- 
sure, P Q . The calculations assume a propellant gas mixture 
of stoichiometric H 2 -O 2 (P c j = latm, T c j = 300 K). In 
each case, the final blowdown pressure at the end- wall in 
the PDRE simulations, P c j , is equal to P a . Note that the 
SSRE performance calculations employ the same constant- 
y assumption used in the PDRE calculations. There are 
four performance curves in the figure: a SSRE with a sonic 
nozzle, a SSRE equipped with an optimized C-D nozzle at 
each pressure ratio, a simple detonation tube (effectively 
a PDRE without a nozzle), and a PDRE fitted with a dy- 
namically optimized C-D nozzle throughout the blowdown 
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Fig. 1 Performance comparison of a pulse detonation rocket 
engine (PDRE) with a conventional steady-state rocket engine 
(SSRE) equipped with both sonic and optimized supersonic 
converging-diverging nozzles. The specific heat ratio, y, is 
held constant in all models. The final blowdown pressure at 
the end-wall in the PDRE is equal to the ambient pressure 
(P C J = Pa)- Propellant mixture: stoichiometric Hj-Oj. Pro- 
pellant initial conditions: P c j = latm, T Cy i = 300K. 

process. 

In examining the figure, a PDRE without a nozzle is per- 
haps best initially compared with a SSRE fitted with a sonic 
nozzle. Note that in this case the simple detonation tube 
provides superior specific impulse, compared to a SSRE 
with a sonic nozzle, over the entire range of pressure ra- 
tios. Note, however, that this advantage is relatively minor 
when expanding to the lowest ambient pressures. More- 
over, in general a PDRE without a nozzle does not compare 
favorably to a SSRE fitted with an optimized C-D noz- 
zle at blowdown pressure ratios above ~ 7. This result is 
primarily due to the largely-choked outflow from a sim- 
ple detonation tube, as there is relatively little conversion 
of thermal energy to directed kinetic energy. However, it 
is also possible to calculate the performance of a detona- 
tion tube fitted with a dynamically optimized C-D nozzle 
throughout the blowdown process. This analysis effectively 
assumes quasi-steady flow within the nozzle at each time 
step in the simulation. The calculations show that if such a 
nozzle could be fitted to a PDRE, then the specific impulse 
of the device would exceed that of a comparable SSRE over 
the entire range of pressure ratios. While such a nozzle is 
a considerable idealization, these simple constant- y results 
clearly demonstrate the critical role that nozzles play in the 
performance of PDREs, particularly at the high pressure 
ratios typical for boost or in-space rocket applications. 

The purpose of this paper is to study the effect of straight 
extensions and C-D nozzles on PDRE gasdynamics and 
performance, in greater detail than can be obtained us- 
ing the previous constant- y method of characteristics ap- 
proach. Details of the quasi- 1-D, finite-rate chemistry com- 
putational fluid dynamics (CFD) model developed by the 
author are provided first. Four different PDRE geome- 
tries are evaluated: a baseline detonation tube, a detonation 


tube with a straight extension, and a detonation tube with 
two types of converging-diverging nozzles. The effect of 
extension length, and C-D nozzle area ratio, on the gas- 
dynamics, performance, and blowdown time of a PDRE is 
studied over a wide range of blowdown pressure ratios ( 1- 
1000). The overall effectiveness of straight extensions and 
converging diverging nozzles in improving performance is 
evaluated. The results are also compared to a steady-state 
rocket system using similar modeling assumptions. 


Governing Equations 

The quasi-one-dimensional Euler equations are the gov- 
erning model used in this study. A quasi- 1-D model is 
at best an approximation of a pulse detonation or steady- 
state rocket engine. Real detonation waves exhibit clear 
multi-dimensional behavior. Additionally, flow separation 
in rocket engine nozzles is strongly dependent on multi- 
dimensional and viscous effects . However, as the primary 
goal of this study is to determine the effect of nozzles on 
pulse detonation rocket systems to first order, the quasi- 1- 
D Euler equations were deemed a reasonable model. The 
time-dependent form of these equations, including the ef- 
fect of finite-rate chemistry, is written as follows: 


au 1 3FS 
dt + s dx 


= H + W 


( 1 ) 


where S represents the cross sectional area of the duct, 
and is a function of distance, x. The state vector, U, and 
convective flux vector, F, are given by 
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The source vectors accounting for the effect of area change, 
H, and for the effect of finite-rate chemistry, W, are 
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Here, pi, . . . , pns represent the mass densities of the ns 
individual chemical species in the problem. The total gas 
density p is sum of the individual species densities, p = 
p*. u is the gas velocity in the x direction, p is the 
gas pressure determined by the ideal gas law. 



where Af* is the molecular weight of the species, R is the 
universal gas constant, and T is the gas temperature, e is 
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the total energy per unit mass of the gas, e = e + ^u 2 , 
where € is the specific energy per unit mass, and for the 
thermally perfect gas model considered here e is a function 
of temperature and gas composition, 

ns 

«-£lta(D (5) 

*= I 

where K* represents the mass fraction of species k, Yt = 
Pk/P , and (k(T) is the specific energy per unit mass of 
species k, as a function of temperature. In this work, the 
ideal gas thermodynamic fits of McBride et al. 12 are used 
to calculate the energy, specific heat, and Gibbs free energy 
for the 9 species in the problem (N 2 , 02 , H 2 , OH, H 2 O, H, 
O, HO 2 , and H 2 O 2 ). 


Chemistry Model 

The source terms for finite-rate chemistry, in Eq. 3, 
are calculated from the sum of relevant chemical reactions 
for each species. Thus, for a chemical mechanism of nr 
elementary reactions, each reaction is expressed as 

ns ns 

(6) 

/=! /=! 


where the v' l n are the stoichiometric coefficients of species 
/ on the reactant side of reaction n, while the v" n are the 
corresponding stoichiometric coefficients on the product 
side. The production term for each species may be writ- 
ten as 


w * s £(<» - V U 

n=l 




and k / and kb represent the forward and backward reaction 
rates, respectively. Details of the chemical kinetics mecha- 
nism for H 2 -O 2 combustion used in this work are shown in 
Table 1. This mechanism was developed by Petersen and 
Hanson 13 to model H 2 -O 2 ignition at the elevated pressures 
typical of practical high-speed propulsion systems. The 
rate coefficients for reactions without pressure dependence 
take on the conventional Arrhenius form: 


k f (T) = AT n cxf>(-E M /RT) (8) 

where A is a constant, and E^t is the activation energy in 
cal/mol. All reactions in the mechanism are reversible. The 
reverse reaction rate, kb, is calculated from the forward rate 
(Eq. 8) and the equilibrium constant. 

Pressure-dependent reactions are modeled using the ap- 
proach given by Kee: 14 


*/=*oo[/V/(l + / > r)]F 
where the reduced pressure, P r , is 
„ *o[M] 


(9) 

GO) 


and the correction factor, F, is in the Troe 15 form: 


InF 


= ( 


1 + 


In P r +c 


1-1 


In F c 


n —d(lnP r + c) 

The Troe centering parameter, F c is given by 
F c = (1 — a)exp(— r/r***) + aexp(— r/r*) + 


(11) 


exp(— r*/7’) 


( 12 ) 


where the constants c = —0.4 — 0.67 In F c , it = -0.75 — 
1.27 In Fc, and d = 0.14. 


Numerical Method 


The governing quasi- 1-D Euler equations are solved in 
finite-volume form throughout the entire domain. The 
equations are solved in split form, in which a fluid convec- 
tion subroutine solves Eq. 1 without the chemistry source 
vector W, 


aU 1 9FS 

aT + s 17 =h 


(13) 


and a finite-rate chemistry integration subroutine solves 
Eq. 1 as a system of ordinary differential equations ignor- 
ing F and H, 



(14) 


The code utilizes the Strang 16 second-order time step 
splitting approach to couple the fluid convection and finite- 
rate chemistry subroutines. For each complete time step, 
the chemistry subroutine is first called for one-half of a 
time step, followed by a full fluid convection time step, and 
then followed by another half-step of the chemistry routine. 
Thus, in operator notation, the solution U at time-step n + 1 
is 

U n+1 = Cc ,/2 CfCc ,/2 U" (15) 

where £/ represents the fluid subroutine solving Eq. 13 
and £ c represents the finite-rate chemistry subroutine solv- 
ing Eq. 14. The time step splitting approach allows both 
high-quality fluid and chemical solvers to be developed and 
tested independently, and then joined together in relatively 
straightforward fashion. As described by Oran and Boris, 17 
time-step splitting works well when relatively small time- 
steps are used. Since the global time-step in the present 
code is strictly limited by the explicit Courant-Friedrichs- 
Lewy (CFL) condition, this is not thought to be a problem. 

The fluid solver used here is the explicit, 2nd-order ac- 
curate (in time and space), symmetric-TVD algorithm de- 
scribed by Yee. 18 The solver employs Roe’s approximate 
Riemann solver modified for nonequilibrium ideal gases. 19 
It also incorporates suggestions by Larrouturou 20 to ensure 
species positivity. 

As with most reactive flow problems, the time- 
integration of the chemistry mechanism in this work re- 
quires a stiff-ODE solver to ensure accuracy. The method 
used here employs a semi-implicit trapezoidal scheme as a 
predictor. 



- U") = A/W" 


(16) 
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Table 1 Chemical kinetics mechanism for H2-O2 combustion, from Petersen and Hanson. 13 


Number 

Reaction* 


A b 

n b 

E b 

1 

o + h 2 ^h + oh 


5.00 X 10 4 * 

2.70 

6290 

r 

H + O 2 4- M ^ H0 2 4- M 


2.80 x 10 18 

-0.90 

0 

3 

H + O 2 + O 2 ^ HO 2 + O 2 


3.00 x 10 20 

-1.70 

0 

4 

H + O 2 + H 2 0 ^ HCh + H 2 0 


9.38 x 10 18 

-0.80 

0 

5 

H 4- O 2 + N 2 ^ HO 2 4* N 2 


2.60 x 10 19 

-1.20 

0 

6 

H + (>2^0 + OH 


8.30 x 10' 3 

0.00 

14413 

7 

H 4* HO 2 ^ O 2 4- H 2 


2.80 x 10 13 

0.00 

1068 

8 

H4-H02^0H + 0H 


1.34 x 10 u 

0.00 

635 

9 

H4-H 2 0 2 ^H0 2 + H 2 


1.21 x 10 7 

2.00 

5200 

10 

oh + h 2 ^h 2 o+h 


2.16 x 10 8 

1.50 

5200 

\\ de 

0H + 0H4-M^H 2 02 + M 

koo 

7.40 x 10 13 

-0.40 

0 



*0 

2.30 x 10 18 

-0.90 

-1700 

12 

0H + H02^02 + H 2 0 


2.90 x 10 13 

0.00 

-500 

13f 

OH-|-H 2 O2^HO 2 4-H 2 0 

ka 

1.75 x 10 12 

0.00 

320 



4 

5.80 x !0 U 

0.00 

9560 

14* 

H0 2 + H0 2 ^0 2 4-H 2 02 

kc 

1.30 x 10“ 

0.00 

-1630 



4 

4.20 x 10 u 

0.00 

12000 

15* 

0 + 04-M^024-M 


1.20 x 10 17 

-1.00 

0 

16* 

04-H + M^OH + M 


5.00 x 10 17 

-1.00 

0 

17 * 

H4-0H + M^H 2 04-M 


2.20 x 10 22 

-2.00 

0 

18> 

H4-H + M^H 2 + M 


1.00 x 10* 8 

-1.00 

0 


Note: species are N 2 , Oj, H 2 , H 2 0, OH, H, O, HO 2 , H 2 O 2 . 
a A11 reactions are reversible. 

b k/(T ) = AT n exp(— E^t/RT); units are in cal, mol, cm 3 , and s. 
C M does not include O 2 , H 2 O, or N 2 ; all collision efficiencies = 1.0. 
Collision efficiencies for M; H 2 = 2.0, H 2 O = 6.0, all others = 1 .0. 
e Troc parameters: a = 0.7346, T*** = 94, T* = 1756, T** =5182. 
f Rate coefficient is non-Arrhenius: k\$ = k a + 4 
g Rate coefficient is non- Arrhenius: *14 = k c +kj 
* Collision efficiencies for M; H 2 = 2.4, H 2 0 = 15.4, all others = 1.0. 

1 Collision efficiencies for M; H 2 = 0.73, H 2 O = 3.65, all others = 1.0. 
^Collision efficiencies for M; H 2 = 1.7, H 2 O - 7.0, all others = 1.0. 


where I is the identity matrix, J is the Jacobian of W, J = 
3 W/dU, and U" +1 is the predicted next value of U. Newton 
iteration of the trapezoidal scheme is used as a corrector, 

[l-±Afj](U? +1 -U£ +, )== 

[U B - U£ +1 + i Af(W B + W£ +1 )] (17) 

where U" +1 is the corrected next value of U. 

The chemistry integration subroutine is designed to use 
several sub-steps, if necessary for accuracy, within a given 
global half-step. The sub-steps are chosen such that no 
species concentration is predicted to change by more than 
1% within each sub-step. Naturally, this leads to more ex- 
tensive use of sub-steps when the concentration within a 
computational cell is far from equilibrium. The Jacobian J 


is evaluated analytically once per sub-step at U n . Thus, the 
LU decomposition of the bracketed matrix terms on the left 
hand side of Eq. 16 can be re-used in Eq. 17. Two corrector 
iterations were performed per sub-step. 

PDRE Simulations 

The four different PDRE geometries studied in this work 
(Fig.2) are highly idealized. Each is based on a constant- 
area (2.0 cm diameter) detonation tube 20 cm in length. 
One end of the detonation tube is closed, while the other is 
either open to the environment (geometry A), attached to a 
straight extension (geometry B) or a converging-diverging 
nozzle section (geometries C and D). Geometry A repre- 
sents the baseline case. Four different straight extensions 
for geometry B were evaluated in this study: L e = 4.0 cm, 
10 cm, 20 cm, and 40 cm. These extensions correspond to 
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ratios of the total assembled tube volume, Vo, to the deto- 
nation tube volume, V t , of Vb/ V, = 1.2, 1.5, 2.0 and 3.0. 
In both of the converging-diverging nozzle geometries (C 
and D), the half-angle in the converging section is 15°. The 
difference between the two geometries is that in C the con- 
verging section is 0.4 cm in length, while in D this length 
is 1.2 cm. This results in a nozzle throat diameter for ge- 
ometry C of 1.786 cm, and for D of 1.357 cm. Thus, the 
maximum rate of outflow is more restricted in geometry D 
than in C, and in both cases is less than that of baseline 
geometry A. The half-angle in the diverging section of the 
nozzle is fixed at 15° in both geometries. The length of 
the diverging section. La, is a variable set at run-time, as 
detailed in the results and discussion to follow. The C-D 
nozzle area ratio is defined throughout this work as the ra- 
tio of the nozzle area at the exit plane to the throat area. 

The detonation tube is pre-filled with a gaseous propel- 
lant mixture with no initial velocity. Stoichiometric H 2 -O 2 
at an initial pressure of 1 atm, and initial temperature of 
300 K is utilized for all calculations in this work. An ideal- 
ized, massless diaphragm isolates the propellant mixture 
from the nozzle and ambient environment until ruptured 
by the detonation wave. The nozzle section is initially 
filled with H 2 gas at the specified ambient pressure. The 
choice of the ambient gas has a nontrivial impact on the 
impulse of geometries B, C and D in single-shot simula- 
tions. Moreover, a number of options can be reasonably 
selected. Hydrogen is used here primarily due to its low 
molecular weight, which allowed a wide range of ambient 
conditions to be studied using the thermochemical model 
explained previously. The impact of using other possible 
choices will be addressed in a limited form in the results 
and discussion section. 

End-wall 

V 

A) 2 cro^P" 



- L t = 20 cm - 


Jdealized 

"diaphragm 


C) 


0.4 cm* + 
— 



Convergences 

half-angle = 15° ^ 

half-angle = 15° 


1.2 cm 


D) 



Convergence 
Ivalf-annla — n- V 


half-angle = 15 Divergence 

half-angle = 15° 


Fig. 2 Schematic of PDRE geometries studied. 


Detonation initiation is facilitated by calculating the 
equilibrium combustion conditions (computed at constant 
p and €) for stoichiometric H 2 -O 2 at 10 atm and 3000 K, 
and then specifying these conditions in the cell adjacent to 
the closed end- wall. 

Ghost cells are utilized to specify the boundary condi- 
tions in the problem. A reflection-type boundary condition 
is utilized at the closed end of the detonation tube to simu- 
late a solid wall. The method of characteristics approach of 
Poinsott and Lele 21 is used to calculate the exit flow bound- 
ary condition. For sonic (choked) or supersonic exit flow, 
all exit flow properties are determined by the interior flow. 
For subsonic exit flow, the ambient pressure is specified 
and the method of characteristics is used to compute the re- 
maining flow properties. If reverse flow is detected at the 
exit, a subsonic inflow boundary condition is calculated us- 
ing the ambient H 2 gas as a constant pressure and enthalpy 
reservoir. 

Two additional special restrictions are imposed in the 
simulations. In order to simulate the effect of an ideal- 
ized diaphragm, only the detonation tube portion of the 
domain (from the closed wall to the diaphragm location) 
is computed initially. A reflection-type boundary condition 
is specified at the diaphragm location until the pressure in 
the adjacent cell rises 1.0% above the initial fill value. This 
special restriction is subsequently removed, and the entire 
domain is computed. Additionally, there is a check per- 
formed when the exit flow is supersonic. Since in this case 
the exit boundary conditions are entirely calculated from 
the interior flow, there is no way for the exit flow to return 
to a subsonic condition. 8 Therefore, at each time step a 
check is made to determine if the pressure from standing 
normal shock at the exit is less than the ambient pressure. 
If true, then the normal shock properties are specified in 
the last interior cell, and a subsonic outflow boundary con- 
dition computed at the exit. 

Grid Resolution 

Grid resolution is a critical issue in all finite-rate chem- 
istry CFD models. Due to the extremely fast kinetics of 
the H 2 -O 2 system, numerical modeling of the reaction zone 
of undiluted H 2 -O 2 Chapman- Jouguet (C- J) detonations re- 
quires an exceptionally fine grid spacing. The C-J detona- 
tion velocity of stoichiometric H 2 -O 2 at an initial condition 
of 1 atm and 300K is calculated by the NASA CEA ther- 
mochemical code 22 to be 2837.1 m/s. Using the calculated 
frozen shock properties at this velocity, and the chemi- 
cal kinetics mechanism described previously, the length of 
the C-J reaction zone for this mixture can be estimated as 
~ 0.02 mm. Thus, a grid spacing of roughly 1 /xm would 
be required to resolve the reaction zone. This is not cur- 
rently practical, nor is it necessary for the performance 
estimates of this study. The chief requirement here is that 
the grid resolution accurately calculate the C-J state and 
velocity. 

A grid resolution comparison of pressure and tempera- 
ture profiles of a propagating C-J detonation, t = 30 /zs 
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Fig. 3 Grid resolution comparison showing pressure (upper 
panel) and temperature (lower panel) profiles of a propagat- 
ing C-J detonation wave at / « 30 /is after initiation. 


after initiation, is shown in Fig. 3. The resolutions vary 
over one order of magnitude. As even the finest resolution 
here, Ax = 0.04 mm, is twice the estimated reaction zone 
size, none of the resolutions can accurately capture the ig- 
nition process behind the leading shock wave. Rather, in 
all cases, the propagating C-J detonation is smeared over 
several cells, resulting in the pressure spike and subsequent 
relaxation toward the C-J state in the finer grids. All four 
resolutions capture the detonation velocity with reasonable 
accuracy (~ 2836 m/s for Ax = 0.04mm, ~ 2840 m/s for 
the other resolutions), though the precision of this estimate 
is naturally lower for the coarser grids. Note, however, that 
in two finer grids the pressure relaxes to the C-J value in 
a significantly shorter distance than in the coarser grids. 
Based on these observations, a uniform grid spacing of Ax 
= 0.10 mm was utilized for all simulations in this study. 

Thrust Calculations 

The time-dependent thrust is calculated at each time step 
by two different methods. One measure of the thrust as- 
sumes a control volume tightly bounding the solid surfaces 
of the PDRE, and is determined by integrating the pressure 
difference across all surfaces. A second measure assumes a 
rectangular control volume encapsulating the PDRE. This 
measure of thrust is determined from the sum of the time 


rate of change of the internal momentum integral across the 
domain, the momentum flux from the nozzle section, and 
the pressure difference across the control volume. In gen- 
eral, there is excellent agreement between the two thrust 
calculations, and the time-integrated impulse calculations 
agree to within 0.1%. In all simulations, the calculation 
proceeds until the pressure at the closed end-wall is equal 
to the ambient pressure ( P c j = P a )- Thus, these simula- 
tions should be thought of as single-shot results. 

Results and Discussion 

There are several critical factors which can impact the 
performance of a PDRE system. The blowdown pressure 
ratio is one of the most important parameters governing the 
performance of any rocket-type system. Additionally, the 
length of a straight extension, or the area ratio of a C-D noz- 
zle, can have a large effect on PDRE performance. In the 
following section, the essential gasdynamics of the blow- 
down process in the four PDRE geometries are reviewed. 
This is followed by a parametric study of the effect of blow- 
down pressure ratio, extension length, and C-D nozzle area 
ratio on specific impulse and blowdown time. The rela- 
tive effectiveness of straight extensions and C-D nozzles 
in improving performance is compared. A comparison of 
the performance of a PDRE fitted with an optimized C-D 
nozzle with an equivalent steady-state rocket engine is also 
provided. 

PDRE Blowdown Gasdynamics 

The blowdown process of a PDRE can be thought of as 
occurring in several different phases, typically differenti- 
ated by the location of important wave structures in the 
tube. Consider the simplest geometry, the baseline deto- 
nation tube (geometry A), at a blowdown pressure ratio of 
10 ( P a = 0.1 atm). Pressure profiles, at several representa- 
tive points in the blowdown history, are shown for this case 
in the upper panel of Fig. 4. The detonation is initiated at 
the end- wall at t = 0. The subsequent propagation of the 
detonation wave is clearly shown in the pressure profile la- 
belled t = 60 /xs. The sharp pressure rise associated with 
the detonation front is followed by a Taylor expansion wave 
which acts to slow and expand the combustion products 
behind the detonation front Well behind the detonation 
front, the products reach a steady pressure at zero for- 
ward velocity. Once the detonation reaches the exit of the 
tube and breaks the idealized diaphragm, outflow begins, 
and an expansion wave propagates back upstream. This 
expansion wave expands and accelerates the combustion 
products toward the exit, and is visible in the t = 180 /xs 
profile. Once this expansion reaches the end-wall, it re- 
flects and begins propagating back downstream toward the 
exit again (f = 300 /xs curve). Eventually, the expansion 
wave reaches the tube exit, and further reduces the pres- 
sure at that location ( t = 420 /xs). Note that the pressure 
at the end- wall is steady at a “plateau” value (P 3 ) until the 
primary expansion wave reaches the end-wall. Similarly, 
the pressure at the exit plane also rapidly achieves a quasi- 
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Fig. 4 Pressure profiles for PDRE geometry A (upper panel) 
and B with L$ ■ 10cm extension Gower panel) at representa- 
tive points in the blowdown history. Blowdown pressure ratio 
is 10 (Pa a 0.1 atm). Propellant mixture: stoichiometric Hj- 
O 2 . Propellant initial conditions: P Cj , m 1 atm, T c>j a 300 K. 

steady value, which holds until the primary expansion wave 
again reaches the exit after reflection from the end-wall. 

The addition of an extension to the detonation tube (ge- 
ometry B) primarily serves to increase the amount of time 
required for some of the waves to traverse the tube. Pres- 
sure profiles at representative times in the blowdown his- 
tory of a detonation tube with a L e = 10 cm extension 
section are shown in the lower panel of Fig. 4. As above, 
the blowdown pressure ratio of this example case is 10. As 
in the baseline case, a primary expansion wave is initiated 
when the detonation reaches the end of the detonation tube. 
Additionally, as there is no combustible gas mixture in the 
extension, the detonation wave changes into a transmitted 
shock wave in the extension. When the transmitted shock 
reaches the exit of the extension section, it initiates a sec- 
ondary expansion wave which also travels upstream toward 
the end-wall. This leads to pressure profiles quite simi- 
lar to the baseline case. However, the overall effect of the 
extension is to slow the rate of pressure relaxation at the 
end-wall. 

A converging-diverging nozzle has a significant effect on 
the pressure distribution of a PDRE. Pressure profiles are 
shown for two C-D example cases in the upper (geome- 




Flg. 5 Pressure profiles for PDRE geometry C with Lg m 
3.6cm (upper panel), and D with Lg =* 2.8cm (lower panel), 
at representative points in the blowdown history. Blowdown 
pressure ratio is 10 (Pa ■ 0.1 atm). Propellant mixture: stoi- 
chiometric H 2 -O 2 . Propellant initial conditions: P e g = 1 atm, 
Tc,i »300K. 

try C with Lg = 3.6 cm) and lower panels (geometry D 
with Lg = 2.8 cm) of Fig. 5. The blowdown pressure ra- 
tio of these cases is again 10. As was the case previously, 
a primary expansion wave, initiated when the detonation 
reaches the end of the detonation tube, is a key factor in 
controlling the pressure relaxation at the end-wall. How- 
ever, it is clear from examining he figures that the C-D 
nozzle acts to expand the combustion products to a much 
lower exit pressure than would otherwise be accomplished 
by the detonation tube alone. The primary difference be- 
tween the two C-D nozzle cases is the significantly smaller 
throat in geometry D, which results in reduced mass out- 
flow compared to geometry C. This flow restriction is sig- 
nificant enough that the primary expansion wave is led by a 
weak shock wave, clearly visible in the lower panel of the 
figure. 

The impulse history of all four of these example cases is 
shown in Fig. 6. The instantaneous thrust can be thought 
of as the slope of these impulse curves. As discussed pre- 
viously, the primary effect of the extension in geometry B 
is to reduce the rate of pressure relaxation at the end- wall. 
This results in a modestly higher thrust in the later portions 
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of the blowdown history, compared to the baseline geom- 
etry A. In contrast, the C-D nozzle has a more significant 
effect on the thrust and impulse. Geometry C has a com- 
paratively higher thrust and resultant impulse throughout 
the entire outflow portion of the blowdown history. Due 
to the more restricted mass outflow, and internal shock 
wave reflections, the instantaneous thrust in geometry D 
is generally lower than that of the other cases. However, al- 
though the overall blowdown time is significantly longer 
than geometry C, the final impulse is nearly equivalent. 
One additional important observation which can be made 
from Fig. 6 is that typically the integrated impulse will 
achieve a value quite close to, or even slightly exceeding, 
the final value well before the end of the simulation. This is 
particularly true of the two C-D nozzle geometries, and is 
observed over the entire blowdown pressure range studied 
in this work. This point will be addressed in greater detail 
in the following sections. 

PDRE Performance Results 

The performance trends just discussed can be seen in 
greater detail by studying the performance of the various 
PDRE geometries over a wide range of blowdown pressure 
ratios. The specific impulse of the baseline detonation tube 
(geometry A), and the detonation tube with several exten- 
sion lengths (geometry B with L e = 4, 10, 20, and 40cm) 
is shown in Fig. 7. As is evident from the figure, the gen- 
eral performance trends of all five configurations are quite 
similar. There is a fairly significant improvement in perfor- 
mance (30-38%) over the blowdown pressure ratio range 
of 1-10, while decidedly less so (4-8%) from 10-100. 
There is very little improvement in performance (~ 1%) 
over the blowdown pressure ratio range of 100-1000. In 
general, there is a definite performance improvement pro- 
vided by the extension section, and this benefit increases 
with extension length. The one exception to this trend 
occurs at a blowdown pressure ratio of two. In this par- 
ticular case, there is no performance benefit for extensions 



Fig. 6 Impulse history for example PDRE geometries. Blow- 
down pressure ratio is 10 (P a = O.latm). Propellant mixture: 
stoichiometric H 2 -O 2 . Propellant initial conditions: P c>l - = 
latm, T fjI * = 300K. 



Fig. 7 Effect of blowdown pressure ratio on specific impulse 
for PDRE geometries A and B. Propellant mixture: stoichio- 
metric H 2 -O 2 . Propellant initial conditions: P c ,j = 1 atm, T c j 
= 300K. 

longer than 10cm. This is due to the fact that the pri- 
mary expansion wave, generated at the interface between 
the detonation tube and the extension section, is alone suf- 
ficient to reduce the pressure at the end- wail to the ambient 
value. In this particular case, this event occurs before the 
secondary expansion wave (generated when the transmitted 
shock reaches the exit of the extension) can reach the end- 
wall. Thus, there is no effect from the longer extensions at 
this pressure ratio. At higher blowdown pressure ratios (5- 
1000) this phenomena does not occur; the pressure decay 
(to ambient) provided by the primary expansion at the end- 
wall is long enough that the secondary expansion wave can 
reach the end-wall in time to affect the pressure relaxation 
rate. 

Both the blowdown pressure ratio, and the nozzle area 
ratio, play critical roles in the performance of the C-D noz- 
zle PDRE geometries. Performance results as a function of 
C-D nozzle area ratio for geometries C and D (upper and 
lower panels, respectively) are shown in Fig. 8. In each 
case, the area ratio was adjusted by varying the length of 
the diverging section of the nozzle (Lj). Results for blow- 
down pressure ratios of 1-1000 are plotted. As is the case 
for conventional rocket-type systems, for a given pressure 
ratio, the performance gain provided by a C-D nozzle will 
generally increase with area ratio up to a certain optimum 
point. After this point performance will decrease with in- 
creasing nozzle area ratio. The dashed curve in each graph 
shows the fixed nozzle area ratio which yields optimum 
specific impulse for each blowdown pressure ratio. It is im- 
portant to recognize that this optimum point is specific to 
the C-D nozzle geometry and blowdown time assumption 
used. 

In general, there is little performance difference between 
the two C-D nozzle geometries. The performance of ge- 
ometry D, at a given area and pressure ratio, tends to be 
slightly lower than geometry C, primarily due to losses 
induced by the reflected shock waves present in the blow- 
down history of D. Note that a C-D nozzle is generally 
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Fig. 8 Effect of C-D nozzle area ratio on specific impulse 
for PDRE geometry C (upper panel) and D (lower panel). 
The dashed line in each graph represents the fixed nozzle 
area ratio giving optimum performance at each pressure ra- 
tio. Propellant mixture: stoichiometric H 2 -O 2 . Propellant 
initial conditions: P C( < = 1 atm, T c>< a 300 K. 


more effective at improving performance at higher pres- 
sure ratios than the extensions in geometry B. The specific 
impulse using an optimized C-D nozzle improves by 53- 
57% over a blowdown pressure ratio range of 1-10, and 
20-21% over 10-100. Even over a blowdown pressure 
ratio range of 100-1000, the performance increases 10%. 
This is largely due to the ability of a C-D nozzle to convert 
the thermal energy of the combustion products into directed 
kinetic energy. This effect can also be seen by plotting spe- 
cific impulse for all four geometries as a function of the 
ratio of total assembled tube volume to detonation tube vol- 
ume, Vo/ V, (Fig. 9). The blowdown pressure ratio is 1000 
in the upper panel of the figure, and 10 in the lower panel. 
It is evident from the figure that a C-D nozzle is generally 
a more effective use of volume than a straight extension, 
particularly at relatively high blowdown pressure ratios. It 
should be noted, however, that this observation is based on 
the assumption that the volume ratio range under consider- 
ation is within the range where increased nozzle expansion 
provides improved performance. 



Fig. 9 Specific impulse plotted as a function of the ratio of to- 
tal assembled tube volume to detonation tube volume, Vj/Vj. 
Blowdown pressure ratio = 1000 (upper panel) and 10 (lower 
panel). Propellant mixture: stoichiometric H 2 -O 2 . Propellant 
initial conditions: P f> / = 1 atm, T e>i a 300 K. 

PDRE Blowdown Time Results 

As the cycle time of an operational PDRE system is of 
great practical interest, it is useful to examine the effect of 
the blowdown pressure ratio, straight extension length, and 
C-D nozzle area ratio on the blowdown time of the PDRE 
geometries considered here. As has been discussed previ- 
ously, the stopping point for all simulations of this study is 
the time at which the pressure at the end- wall decays to the 
ambient value. For the purpose of this discussion, we will 
define the blowdown time, tbd, as equal to this time value. 
The strong dependence of tbd on blowdown pressure ratio 
and extension length can be seen in Fig. 10. The time in 
that figure has been normalized by Uc-j/L, = 70.495/xs, 
which is the characteristic time feu- a C-J detonation to tra- 
verse the detonation tube. Note that the highest blowdown 
pressure ratios, there is an almost directly proportional de- 
pendence of blowdown time on extension length. The ob- 
vious kink in the L e = 20 cm and L e = 40 cm curves at a 
blowdown pressure ratio of two was explained in the previ- 
ous section. 

The blowdown time for the two converging-diverging 
nozzle PDRE geometries is shown in Fig. 11. The val- 
ues have again been normalized by Uc-j/L,. As would 
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be expected, tbd is consistently longer for geometry D than 
C, over the entire range of blowdown pressure ratios. Note 
that these results assume the use of the optimum fixed C- 
D nozzle for each pressure ratio, though since the nozzle 
throat is largely choked, tbd is nearly independent of Lj at 
a given pressure ratio. One of the notable results of these 
simulations is that the instantaneous thrust of even an op- 
timum fixed C-D nozzle is slightly negative for much of 
the latter portions of the blowdown time. This result is due 
to the fact that the C-D nozzle overexpands the residual 
combustion products when the pressure in the detonation 
tube has decayed sufficiently. Thus, it is useful to also 
plot the time at which the thrust first drops to zero in the 
blowdown history (f z/ ). It is evident from the figure that 
this point occurs considerably earlier than tbd* particularly 
at the highest pressure ratios. As the instantaneous thrust, 
even while negative, is so small at this point, there is lit- 
tle difference in integrated impulse between these two time 
points (typically at most a loss of ^ 1 s of specific impulse 
over the interval). However, this observation suggests that 
a considerable savings in overall cycle time could be real- 
ized with an actual slight gain in performance. Note that as 
the end- wall is the only thrust surface for geometries A and 
B, t zt = tbd for those cases. 

Equivalent SSRE Model and Performance Comparison 

It is instructive to compare the specific impulse of a 
PDRE fitted with a C-D nozzle with a SSRE under equiva- 
lent modeling assumptions. While frozen and equilibrium 
rocket performance calculations can be readily obtained 
from the CEA code, these results are not directly compa- 
rable to the finite-rate chemistry model used in the current 
PDRE code. The primary concern is the tendency for 
chemistry to slow down in real nozzle systems as the tem- 
perature and pressure are reduced in the expansion process. 
Thus, it is best to compare the finite-rate PDRE systems 
with a finite-rate SSRE model. This SSRE CFD model 
is heavily derived from the PDRE code. The 20cm det- 



Blowdown Pressure Ratio (P^/ PJ 


Fig. 10 Normalized blowdown time UmUc-j/Li) (dotted as 
a function of blowdown pressure ratio for PDRE geometries A 
and B. Propellant mixture: stoichiometric H 2 -O 2 . Propellant 
initial conditions: P e j = 1 atm, T c> , = 300 K. 



Blowdown Pressure Ratio (P^/ PJ 

Fig. 11 Normalized time ( tUc-j/Lt ) plotted as a function 
of blowdown pressure ratio for PDRE geometries C and D. 
Propellant mixture: stoichiometric H 2 -O 2 . Propellant initial 
conditions: P Ct ; = 1 atm, T c ,j = 300 K. 

onation tube is replaced with a thrust chamber 0.1 cm in 
length, and 4.0 cm in diameter (Fig. 12). This thrust cham- 
ber then converges at a constant 15° angle, over a length of 
4.13 cm, to a throat 1.787 cm in diameter. After the throat, 
the nozzle diverges at 13° to a length specified by the user. 
Thus, the throat diameter and nozzle expansion rate of the 
SSRE model closely matches that of PDRE geometry C. 
The equilibrium temperature, pressure and composition of 
stoichiometric H 2 -O 2 , burned at constant pressure and en- 
thalpy from initial conditions of 300K and 1 atm, are fed 
as a constant enthalpy reservoir inflow boundary condition 
into the domain. Thus, this analysis assumes that the pro- 
pellant feed system of the SSRE model is the same as the 
PDRE. 

Similar to the PDRE nozzle optimization study, the 
finite-rate SSRE CFD model is run at various expansion 
ratios until the optimum specific impulse is obtained for a 
given pressure ratio. However, as the solution converges 
to steady-state in this model, flowfield information from 
the highest pressure ratio case can be used to guide area 
selection at the lower pressure ratios. In each case, the 
SSRE model is run until the solution converges. In general, 
the specific impulse and optimum expansion ratio using 
finite-rate chemistry are slightly larger than that for frozen 
chemistry. Additionally, if the chemistry is frozen in the 
SSRE CFD model, there is excellent agreement (to within 
0.3% in specific impulse) with the frozen-chemistry predic- 
tions of CEA. 

The mixture-based specific impulse for the baseline det- 
onation tube (PDRE geometry A), PDRE geometry C fitted 
with an optimized fixed C-D nozzle (at each blowdown 
pressure ratio), and the equivalent optimized SSRE sys- 
tem is plotted for blowdown pressure ratios of 1-1000 in 
Fig. 13. As would be expected from the results discussed 
previously, the relative gain from the C-D nozzle systems 
becomes more pronounced at higher pressure ratios. It is 
interesting to note that that both the baseline and C-D noz- 
zle PDRE systems outperform a SSRE at pressure ratios 
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below ~ 7. This is due to the fact that the relative gain 
from the extra pressurization provided by the detonation 
process is greatest in this regime. Thus, the PDRE may 
have considerable theoretical potential for rocket-type ap- 
plications when the pressure of the ambient environment 
is high. Additionally, at higher blowdown pressure ratios, 
a PDRE with an optimized fixed C-D nozzle still has a 
greater specific impulse than an equivalent SSRE. These re- 
sults indicate that a PDRE equipped with a fixed C-D noz- 
zle, though obviously not as efficient at generating thrust as 
a (theoretical) dynamically adaptive nozzle, can still yield 
performance superior to an equivalent SSRE over a wide 
range of pressure ratios. Note, however, that this perfor- 
mance gain becomes relatively smaller at higher pressure 
ratios. 


Conclusions 

A quasi- 1-D, finite-rate chemistry CFD model for study- 
ing PDRE gasdynamics and performance is described and 
implemented. The single-shot performance and blowdown 
time characteristics of four different PDRE geometries are 
studied, over a range of blowdown pressure ratios from 1- 
1000. The pressure at the end-wall is allowed to decay 
to the ambient value for all simulations in this work. The 
results show that both straight extensions and converging- 
diverging nozzles can provide valuable improvements to 
the performance of a baseline detonation tube. However, 
optimized C-D nozzles are generally more effective in in- 
creasing performance than straight extensions, particularly 
at higher pressure ratios. Studies of the blowdown time 
of the PDRE systems show a strong dependence on pres- 
sure ratio for all geometries. Straight extension length will 
also generally increase blowdown time. Examination of 
the blowdown history of the optimized C-D nozzle geome- 
tries shows that the blowdown of these systems could be 
cut off well before the end- wall pressure reaches the ambi- 
ent value, with an actual slight gain in performance. The 
performance of a PDRE fitted with an optimized fixed C-D 
nozzle is also compared with a SSRE using similar mod- 
eling assumptions. The results show that while a PDRE 
enjoys the largest relative performance gain over a SSRE at 
low pressure ratios, there is still a noticeable relative advan- 
tage for the PDRE at a blowdown pressure ratios of 1000. 
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Fig. 12 Schematic of SSRE geometry studied. 



Fig. 13 Performance comparison of various pulse detona- 
tion and steady-state rocket devices. All results obtained using 
finite-rate quasi-l-D CFD calculations. The final blowdown 
pressure at the end-wall in the PDRE is equal to the ambient 
pressure (P Cf f =* Pm)- Propellant mixture: stoichiometric 
H 2 -O 2 . Propellant initial conditions: P Cy i » latm, T Cyi = 
300 K. 
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